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Abstract. We present 7-mm continuum observations of 14 low-mass pre-main-sequence stars in the Taurus- Auriga 
star-forming region obtained with the Very Large Array with ~l'.'5 resolution and ~0.3 mjy rms sensitivity. 
For 10 objects, the circumstellar emission has been spatially resolved. The large outer disk radii derived suggest 
that the emission at this wavelength is mostly optically thin. The millimetre spectral energy distributions are 
characterised by spectral indices a mm = 2.3 to 3.2. After accounting for contributions from free-free emission and 
corrections for optical depth, we determine dust opacity indices f3 in the range 0.5 to 1.6, which suggest that 
millimetre-sized dust aggregates are present in the circumstellar disks. Four of the sources with j3 > 1 may be 
consistent with submicron-sized dust as found in the interstellar medium. Our findings indicate that dust grain 
growth to millimetre-sized particles is completed within less than 1 Myr for the majority of circumstellar disks. 
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1. Introduction 

Although dust grains constitute only a minor fraction 
of the disk material around young stars (~1 Myr), they 
play a pivotal role in the complicated multi-stage pro- 
cess of planet formation. The growth of submicron-sized 
dust particles as found in the interstellar medium (ISM) 
to larger aggregates by coagulation (collisional sticking) is 
the first essential step towards building larger rocky bod- 
ies (planetesimals) that may event ually accrete into terre s- 
trial planets and planetary cores ijBeckwith et alJ feoOOI. 

Statistical analysis of the near-infrared excess of 
young stellar clust ers shows that disk emission disappear s 
within a few Myr l|Strom et allll989t lHaisch et aljhoOljh 
Similarly, sub-mm and mm observations suggest a de- 
crease in the amount of cold circumstellar d ust during the 
post T- Tauri phase l|Carpenter et al J 12*0051. The decline 
of dust emission can be interpreted as a consequence of 
gradual mass loss through disk dissipation and/or opacity 
changes due to particle growth. It is important to stress 
that such studies only trace the evolution of the dust disk. 
The evolution of the gaseous disk, containing the bulk of 
the disk mass, is at present only poorly understood. 

For many years there has been little sound observa- 
tional evidence for dust grain growth in disks around pre- 
main-sequence stars. The analysis of spectral energy dis- 
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tributions showed that the sub-mm/mm fluxes of T Tauri 
stars decline more slowly towards lon ger wavelengths than 
expec ted for ISM-sized dust (e. g. iBeckwith fc SargenlJ 
1991). It was tempting to interpret a shallow spectral 
slope as an indication of the presence of particles much 
larger than in the ISM (^> 0.1 /im). It was soon realised, 
however, that spatially unresolved disk observations can- 
not be used to distinguish between small, optically thick 
disks containing sub-micron ISM dust and extended, op- 
tically thin disks with larger dust particles. Spatially re- 
solved jmaj!£S_arej^eded to^ degen- 
eracy l|Koerner et alJll995t IPutrev et alJll99fj|) . 

Recent technical improvements at the Very Large 
Array (VLA) allowed the resolution of the disk around 
the young star TW Hya at 7 mm and the determination 
of the dust opacity index jWilner et al l2000HCalvet et alJ 
120021) . iTesti et alJ l)2003jl resolved the dusty disk around 
the pre-main-sequence star CQ Tau and concluded that 
millimetre- and even centimctrc-sizcd dust particles must 
be present in the outer disk. In a similar stu dy of six iso- 
lated intermediate- mass (Herbig Ae) stars bv lNatta et alJ 
(2004), dust grain growth was inferred for two objects. 

In this paper, we report the results of 7-mm continuum 
observations of 14 T Tauri stars located in the Taurus- 
Auriga star-forming region, currently the largest sample of 
low-mass pre-main-sequence stars investigated for signs of 
dust grain growth to millimetre/centimetre particle sizes. 
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Table 1. Properties of sample stars. 



Star 


Coordinates (J2000.0) 


Spectral 


T ( „ b) 

cff 


Mi c) 


ii b) 


Age (c) 


Location' d ' 




„ /h m s\ 
V / 


o ( ) 


type (a) 


(K) 




v (jy 


(Myr) 




CY Tau 


4 17 33.74 


+28 20 46 5 


Ml 


3720 


0.48 


0.47 


0.7 


L1495 


RY Tau 


4 21 57.42 


+28 26 35 5 


Kl 


5080 


1.69 


7.60 


0.2 


L1495 


FT Tau 


4 23 3Q 1 Q 


+24 56 14 9 


cont. 


3890 (c) 


— 


— 


— 


B217 


DG Tau B 


4 27 09 56 

'-t Zi I UZi . iJU 


+96 05 30 4 
n^^u uu ou.4 


— 


— 


— 


— 


— 


B217 


DG Tau 


A 07 04 fiQ 


_i_9fi Dfi 1fi 1 


M 


3890 (c) 


0.56 




0.3 


B217 


HL Tau 


4 31 3£ A 1 


-j-lo lo f .0 


K7-M2 


4060 


0.55 


0.9 (c) 


1.0 


L1551 


GG Tau 


A 39 30 30. 


j_i 7 3i /in i 


K7 


4060 


0.65 


1.50 


0.3 


B217 


UZ Tau E 


4 32 43.10 


+25 52 31.0 


Ml-3 


3680 


0.44 


2.94 (c) 


<0.1 


B217 


DL Tau 


4 33 39.09 


+25 20 37.9 


K7 


4060 


0.56 


0.77 (c) 


1.2 


TMC2 


DM Tau 


4 33 48.74 


+18 10 10.0 


Ml 


3720 


0.62 


0.25 


3.2 


L1551 


CI Tau 


4 33 52.03 


+22 50 30.2 


K7 


4060 


0.70 


0.87 


0.8 


L1536 


DO Tau 


4 38 28.59 


+26 10 49.4 


MO 


3850 


0.72 


1.20 


0.6 


TMC1 


LkCa 15 


4 39 17.81 


+22 21 03.6 


K5 


4350 




0.74 


1.0 (c) 


L1536 


GM Aur 


4 55 10.98 


+30 21 59.5 


K3 


4730 


0.72 


0.83 


1.8 


L1517 



Re ferences — (a ) jHerbig fc BelJ lll988ft : fbl lKenvon fc Hartmannl dl995l) : (c) lBeckwith et alJ Jl99(t) : fdl lLeinert et alJ l|l99" 
fel lDuvert et alJ (|200(t . 



2. Observations 

2.1. Sample of objects 

Our sample consists of 14 low-mass pre-main-sequence 
stars located towards the Taurus- Auriga star formation 
region (Table All sample stars are catalogued as clas- 
sical T Tauri stars, based on the equivalent wi dth of their 
Ha emission line £10 A l|Herbig fc Belli 1988|) . They are 
classified as Class II young stellar objects, as derived from 
the slope of their infrared spectral energy distributions, 
otm = dlogfAFxl/dlogA, with -2 < a < lAdams et all 



IS iKenvou fc Ha.rtma.uulll9Ql iHarlmanrto):* 

Observations of the millimetre continuum emission 
have confirmed the presence of circumstellar disks around 
these T Tauri st ars, with disk masses ranging from 
0.02 to 0.7 M (|Beckwith et all ll99(A iDutrev et aD 
Il996|) . For the majority of the stars in our sample, 
emission from circumstellar CO gas has been detected, 



oitcn tracing a cusk m tvepierian rotation (rvocrncr 
19931 IDutrev et all Il994l iKoerner & Sargentl 


et^j. 
1995: 


Handa et al.lll995: Dutrev et all 


1996; 


Jensen et alJ 


1996; 


Mitchell et all Il997l iDuvert et all 1200(1 iNaiita et all 



2003) 

The stellar ages and masses have been determined 
by fitting pre-main-sequence isochrones to loci of the 
stars in the Hertzsprung-Russell diagram; they range 
from 0.1 to 3.2 Myr and from 0.4 to 1.7 M q, respectively 
l|Beckwith et al.ll990HDuvert et alJl2000|) . The uncertain- 
ties of the age determinations arc considerable and pre- 
clude a meaningful correlation with the opacity indices (3. 

The average distance to the Taurus-Auriga molecu- 
lar cloud as found in th e literature va ries between 135 p c 
i Cernicharo et al.lll985h and 160 pc llStrom et alJll989h 
We adopt the standard value of 140 dc l(Eliasl Il978l 



iKenvon et aDll994l: lLoinard et aTll2005|) to convert angu- 
lar scales into physical sizes. 

2.2. 7-mm continuum observations 

The NRAO 1 Very Large Array was used to observe the 
dust continuum emission of 14 T Tauri stars at 43.34 GHz 
(7 mm). The total receiver bandwidth was 100 MHz, di- 
vided into two 50 MHz bands centred at 43.3149 GHz and 
43.3649 GHz. The observations were carried out in D con- 
figuration on 22-23 March 2003 and again on 12 May 2003. 
The 27 VLA antennas provided 27 x (27-l)/2 = 351 
baselines from 35 m to 1.03 km, corresponding to 5-150 kA 
at 7 mm. The May observations took place during array 
reconfiguration (D to A); only 26 antennas were available. 

Phase calibration was accomplished by sandwich- 
ing the on-source observations between pointings of 
nearby secondary calibrators (radio sources J0431+2037, 
J0403+2600, J0443+3441). The typical source/phase cal- 
ibrator cycle time was 5-10 minutes. Absolute flux cal- 
ibration was obtained from observation of the quasar 
J0542+4951 (3C147) which is assumed to have a flux 
density of 0.91 Jy. The estimated uncertainty in the ab- 
solute flux calibration is ^10%. Reference pointing mea- 
surements were made approximately every hour. 

The (u, v) data sets for the three observing dates were 
merged using the ALVS task DBCDN. CLEANed images 
were obtained with IMAGR using a ROBUST parameter of 
0, which is intermediate between natural and uniform 
weighting of the visibility data and optimises for spatial 
resolution and sensitivity. Three weak sources (DM Tau, 
CI Tau, and LkCa 15) were CLEANed using natural 
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weighting (R0BUST=5) to increase signal-to-noise ratio and 
allow secure detection. The extended emission of GG Tau 
is also best imaged using natural weighting. Figure^Jshows 
the corresponding contour plots for the 7-mm images. 
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Fig. 1. VLA D-configuration images of the A = 7 mm continuum emission. Contour levels are drawn in steps of the 
corresponding 2a rms noise level; the average la rms noise level is ^0.15 mJy/beam. The synthesized beam is shown 
in the bottom left corner of each panel; the average beam size is ~l'/5. 
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Fig. 2. VLA D-configuration images at A = 1.3 cm (left column), A = 2.0 cm (middle column), and A = 3.6 cm (right 
column) for RY Tau, DG Tau and DG Tau B, HL Tau, and UZ Tau. Contour levels are given for each panel separately 
in multiples of 2a rms noise level. The average la rms noise levels are approximately 50, 80, and 30 /iJy/beam at 1.3, 
2.0, and 3.6 cm, respectively. About 30" to the east of HL Tau lies XZ Tau (indicated by crosshairs). 
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2.3. Centimetre observations 

Five of the sources have poor-qu ality archival VLA de- 
tections at 8 GHz llWiluerj 1998). All five sources were 
reobserved at 3.6 cm (8.46 GHz), 2.0 cm (14.94 GHz), 
and 1.3 cm (22.46 GHz). The total bandwidth was again 
100 MHz, with two 50 MHz channels 50 MHz apart. The 
1.3-cm observations were carried out on 17 March 2003; 
the 2.0- and 3.6-cm observations on 21 March 2003. All 
observations were made in the compact D configuration, 
covering baselines of 3-77 kA, 2-50 kA, and 1-30 kA in the 
(u,w)-plane at 1.3, 2.0, and 3.6 cm, respectively. 

The flux density scale was set by observations of the 
calibrator J0542+4951 with known flux densities of 1.80, 
2.71, and 4.74 Jy at 1.3, 2.0, and 3.6 cm, respectively. 
CLEANed images were produced with IMAGR and natural 
weighting of the visibilities (R0BUST=5). Finally, primary 
beam corrections were performed using the AXVS task 
LTESS to compensate for the non-uniform single-antenna 
response across the field-of-view of the entire array. The 
contour maps are shown in Figure note that DG Tau 
and DG Tau B appear in the same image. 

2.4. Radio seeing 

Short-term variations in precipitable water vapour lead 
to changes in the path length of an electromag- 
netic wave travelling through the troposphere, causing 
phase fluctuations when observed by radio interferom- 
eters. Tropospheric phase noise negatively affects spa- 
tial res olution and cohere nce, especially at mm wave- 
lengths l|Carilli et alJll99fli> . 

In order to estimate the amount of smearing intro- 
duced by tropospheric phase noise, we compared self- 
calibrated and standard-calibrated images of two bright 
radio sources. Images of self-calibrated sources are cor- 
rected for the radio seeing on short timescales (3.3 seconds 
integration time), and deconvolution will give the intrin- 
sic source size. The same procedure applied to images ob- 
tained from ordinary calibration gives the 'seeing disk'. 
By using two test sources, one unresolved (J0426+2327), 
the other extended (J0412+2305), we checked the flux re- 
duction by phase-noise induced smearing for point sources 
and extended sources, respectively. 

The analysis of the merged 7-mm data showed that the 
overall radio seeing was ~0'.'5. Inspection of the single- 
day data revealed that the seeing on the second day 
(23 March 2003) was the worst (~0'/7), while on 22 March 
and 12 May better water vapour conditions prevailed 
(about Of.'l and ff.'A, respectively). The integrated 7-mm 
fluxes are reduced by 4-5%, and are within 2% for each 
day. 

The seeing at 1.3 cm was about 0'.'5. The peak fluxes 
might be reduced by up to 5%; the integrated fluxes were 
affected less than 1%. Tropospheric seeing effects at longer 
centimetre wavelengths are expected to be negligible. 



3. Results 

3.1. 7-mm fluxes and disk sizes 

We detected the continuum emission from cold circum- 
stellar dust at 7 mm for all 14 T Tauri stars comprising 
our sample. Two sources (DM Tau and LkCa 15) have 
signal-to-noise ratios < 5, and we do not consider these 
sources further in our analysis. The integrated fluxes were 
measured on the combined maps using the AIVS routine 
IMSTAT, and corrected for seeing reduction. In the light 
of the radio seeing analysis, we decided to use the ex- 
cellent conditions of the first-day data set to determine 
source sizes. The task IMFIT was used to check if the 
source was spatially resolved or unresolved by fitting a 
two-dimensional Gaussian function to the emission peak. 
The deconvolved disk sizes are of the order 100-200 AU, 
representing typical dimensions for protoplanetary disks. 
The major and minor axes can, in principle, be used to 
derive the inclination angle of the disks. Their uncertain- 
ties, however, are very large, and the disk orientations are 
poorly constrained. Table|2]summarises the results derived 
from the 7-mm maps. 

The extended dust emission of three sources is of 
particular interest. The map of HL Tau shows a promi- 
nent "finger" pointing in south-western direction. The 
measured position angle of 235° ± 1 0° is roughly co- 
aligned with the cou nter-jet of HL Tau l|Mundt et al . 1988; 
IWilner fc Lall2000lL 

We detected extended dust emission and re- 
solved circumstellar material around the northern bi- 
nary (~( y/3 separation) o f the quadruple system 
GG Tau l(Leinert et al.lll993h . Previous 13 CO line emis- 
sion and 1.4- mm continuum observations demonstrated 
that the circumbinary disk consists of two structural com- 
ponents: an 80- AU wide ring at a radius of ~220 AU, 
and a fainter, more e xtended disk ijDutrev et al.l Il99l 
iGuilloteau et al.lfl999l) . Given the complicated source ge- 
ometry of GG Tau we have not attempted to deconvolve 
a simple Gaussian model from the 7-mm map. 

The quadruple system UZ Tau could be resolved into 
the bright eastern and the (marginally detected) western 
components, both of which are known to harbour a binary 
system. Speckle and direct imaging observations resolved 
UZ Tau W into a binary system wi t h a projected sepa - 
ration of ~50 AU ijOhez et al.lll993l ISimon et al1ll995h . 
Radial-velocity measurements showed UZ Tau E to be 
a spectrosco pic binary with a projected semimaior axis 
of ~0.1 AU llMathieu et al.lll996t IPrato et allEfflj . We 
confirm the findings of IJenseTietaiT <)l996j) . who noted 
a substantial reduction of millimetre emission around 
UZ Tau W compared to the tight binary UZ Tau E. In 
the former system, the two components have a separa- 
tion comparable to a typical protoplanetary disk, leading 
to disk truncation and therefore smaller and less massive 
disks. The two components of UZ Tau E, on the other 
hand, apparently do not a ffect the common cir cumbinary 
disk in which they reside l)Mathieu et alJEoOfl . 
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Table 2. Seeing-corrected 7-mm flux densities, spatial resolutions, and estimated source sizes. The stated uncertainties 
for the flux densities do not contain the uncertainty of the absolute flux calibration. 



Source Integrated flux Disk resolved? Deconvolved source size P.A. 

(mJy) (AU) H 



CY Tau 


1.19 


± 


0.26 


no 






— 










RY Tau 


2.97 


± 


0.29 


yes 


180 


(±30) 


X 


40 (±50) 


140 


± 


10 


FT Tau 


1.62 


± 


0.27 


partially 


170 


(±50) 


X 


— 


60 


± 


20 


DG Tau B 


3.57 


± 


0.29 


yes 


130 


(±30) 


X 


40 (±50) 


50 


± 


30 


DG Tau 


5.30 


± 


0.37 


partially 


110 


(±30) 


X 


— 


10 


± 


10 


HL Tau 


4.92 


± 


0.43 


yes 


140 


(±40) 


X 


80 (±40) 


170 


± 


30 


GG Tau 


3.24 


± 


0.42 


















UZ Tau E 


1.84 


± 


0.28 


partially 


120 


(±50) 


X 




80 


± 


20 


DL Tau 


1.49 


± 


0.25 


yes 


210 


(±60) 


X 


140 (±100) 


110 


± 


40 


DM Tau 


0.74 


± 


0.20 


















CI Tau 


0.76 


± 


0.17 


partially 


220 


(±60) 


X 




40 


± 


10 


DO Tau 


2.47 


± 


0.29 


partially 


60 


(±50) 


X 




70 


± 


30 


LkCa 15 


0.44 


± 


0.17 


















GM Aur 


1.05 


± 


0.23 


yes 


200 


(±50) 


X 


60 (±70) 


10 


± 


20 



3.2. Centimetre fluxes 

Centimetre emission could be detected for all five sources 
observed at 1.3, 2.0, and 3.6 cm (Table OJ. The integrated 
fluxes of UZ Tau at 2.0 and 3.6 cm are below the 4-er 
level. The counter-jet of HL Tau mentioned in Section l3~Tl 
can also be seen in the corresponding 1.3-cm image. Its 
position angle has been measured to 225° ± 20° , in agree- 
ment with the value measured from the 7-mm map. The 
UZ Tau multiple system is not resolved into its eastern 
and western components. 

Radio emission from low- mass pre-mai n-sequence stars 
is a common phenomenon (Giidcl 2002). VLA observa- 
tions of classical T Tauri stars at radio wavelengths re- 
vealed rising spectral indices and lar ge angular sizes, in - 
terpreted as signs of wind emission IjCohen et al.l ll982). 
Observational evidence and theoretical arguments favour 
focussed anisotropic outflows and collimated jets rather 
than uniform , isotropic mass flows f rom T Tauri stars 
l|Cohenlll982l ICohen fc Bieginglll986h . Wind emission is 
likely to contribute to the 7-mm continuum emission and 
has to be subtracted from the measured 7-mm flux den- 
sities in order to accurately derive dust opacity indices 
from spectral slopes at millimetre wavelengths. 



Table 3. Summary of centimetre observations. 1.3-cm val- 
ues corrected for seeing-related flux reduction. The uncer- 
tainty in the absolute flux density calibration is not in- 
cluded. 



Wavelength 


Source 


Integrated 


Detection 






flux (mJy) 




1.3 cm 


RY Tau 


0.92 ± 0.08 


12 




DG Tau B 


1.23 ± 0.09 


14 




DG Tau 


2.17 ± 0.11 


20 




HL Tau 


1.63 ± 0.11 


15 




UZ Tau 


0.77 ± 0.11 


7 


2.0 cm 


RY Tau 


0.63 ± 0.13 


5 




DG Tau B 


0.80 ± 0.13 


6 




DG Tau 


1.33 ± 0.15 


9 




HL Tau 


0.88 ± 0.14 


6 




UZ Tau 


0.48 ± 0.14 


3 


3.6 cm 


RY Tau 


0.31 ± 0.04 


8 




DG Tau B 


0.46 ± 0.05 


9 




DG Tau 


1.27 ± 0.05 


25 




HL Tau 


0.57 ± 0.05 


11 




UZ Tau 


0.10 ± 0.03 


3 



3.3. Free- free emission 

At centimetre wavelengths the radio continuum of young 
stars is expected to be dominated by free-free radia- 
tion (bremsstrahlung) from ionized winds or outflows. 
There are, however, some observational indications that 
dust continuum emission may be a substantial source 
contributing to the centimetre fluxes of T Ta uri stars 
I Rodriguez et ^Il994l IWiluer et al.lll 9961 Eool. 



With the exception of RY Tau, all other sources ob- 
served at 1.3 cm appear to be resolved, at least in one spa- 
tial direction. The position angles roughly match (within 
20°) those derived from the 7-mm images. This finding 
supports the hypothesis that dust emission rather than an 
ionised wind is the main emission mechanism. For HL Tau 
the position angles differ by ^40°, probably as a result of 
the bright counter-jet. 
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Fig. 3. Free-free emission contribution to 7-mm flux densities. The dashed line represents the fit to the centimetre 
fluxes, assuming that 50% of the 1.3-cm emission arises from free-free radiation. The square depicts the corrected 
7-mm flux after subtraction of the estimated free- free contribution (shown as open circle). The thick line shows the 
power-law fit to the millimetre data; the shaded region indicates the uncertainty range for the 1-7 mm slope. Data 
points shortwards of A=7 mm were compiled from the literature (Tabic 0. 



We measured the free-free spectrum for four sources by 
fitting a power law to the centimetre flux values (FigureEJ) ■ 
The measured 2.0 and 3.6-cm fluxes are assumed to orig- 
inate from free-free radiation only. At higher frequencies, 
a reliable quantification of the corresponding contribu- 
tion of thermal dust emission and free-free radiation is 
difficult; we assumed an equal mixture of dust and wind 
emission at 1.3 cm. The spectral index ag (S^ oc v aii ) 
was found to be in the range —0.1 to +0.4. A sym- 
metric, ionized, and opaque wind with constant velocity 
has a spectral po wer in dex of 0.6 llPanagia fc Fellilll975l: 
IWright Rr. FSa,rlowlll975t l01nonlllQ75h . n« = -0.1 is ex- 
pected fo r totally transpare nt (optically thin) free-free 
emission ( Mezger et al.lll967i) . 



For all four objects we extrapolated the free-free ra- 
dio spectrum to A=7 mm to estimate the contribution of 
wind emission to the 7-mm fluxes. Subtracting the esti- 
mated free-free emission from the measured 7-mm value 
one obtains the thermal dust emission. We found that 
about 20% of the 7-mm emission originates from free-free 
radiation; 80% are due to dust continuum emission. For 
sources where no centimetre data were available, we used 
these numbers to correct the 7-mm fluxes. 



3.4. Spectral indices and dust opacity indices 

The outer disk radii derived suggest that the emission 
is optically thin at 7 mm. One can therefore directly 
use the slope of the millimetre spectral energy distribu- 
tion (F v (xv a ) to derive a dust opacity law {n v oc 
llBeckwith fc Sargendll99ll iTesti etal]l2003t iNatta et alJ 
2004|) • Were the millimetre emission completely optically 
thin, i.e. without any contribution from the optically thick 
inner disk, one could simply read off the opacity index 
(3 from the spectral index a via the relation [3 = a — 2. 
There is, however, a non-negligible contribution from the 
optically thick part of the disk to the measured fluxes. 
Knowing this contribution will enable an improved esti- 
mate for the emissivity index (3. 

The revised relation between the opacity (emissiv- 
ity) index [3 and the observed spectral index a is 
given by j3»(a-2)x(l+A), where A is the ra- 
tio of optically thick to optically thin disk emission 
given by A = -y x { (2 - q) ln[(l - p/2) fl}" 1 (Eq. 20 in 
iBeckwith et al.lll99f)|) . Here p, q, and f are the power-law 
exponent of radial disk temperature, the power-law expo- 
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Table 4. Spectral slopes derived from power-law fitting in the range 1-7 mm using literature values (Table [SJ and 
our 7- mm data (Col. 2). For objects without centimetre measurements, a 20% contribution of free- free emission was 
assumed (denoted by brackets in Col. 3). Spectral slopes corrected for free-free radiation are listed in Col. 4. The 
estimated dust opacity indices /3~ (Col. 5) were corrected for optically depth (using A = 0.2), yielding the final 
values (Col. 6). 



Source 


mm slope 


f-f contribution 


f-f corrected 


estimate 


Opacity index 






at 7 mm (%) 


C^corr 


p — Clcorr 2 


/3 = (i + A)/r 


RY Tau 


2.51 ± 0.10 


24 


2.66 ± 0.09 


0.66 ± 0.09 


0.8 ± 0.1 


FT Tau 


2.62 ± 0.17 


(20) 


2.76 ± 0.24 


0.76 ± 0.24 


0.9 ± 0.3 


DG Tau B 


2.88 ± 0.20 


24 


3.17 ± 0.19 


1.17 ± 0.19 


1.4 ± 0.2 


DG Tau 


2.45 ± 0.09 


20 


2.54 ± 0.11 


0.54 ± 0.11 


0.7 ± 0.1 


HL Tau 


2.93 ± 0.10 


23 


3.04 ± 0.12 


1.04 ± 0.12 


1.3 ± 0.1 


UZ Tau E 


2.56 ± 0.08 


(20) 


2.66 ± 0.09 


0.66 ± 0.09 


0.8 ± 0.1 


DL Tau 


2.73 ± 0.11 


(20) 


2.82 ± 0.14 


0.82 ± 0.14 


1.0 ± 0.2 


CI Tau 


3.00 ± 0.27 


(20) 


3.12 ± 0.30 


1.12 ± 0.30 


1.3 ± 0.4 


DO Tau 


2.29 ± 0.06 


(20) 


2.38 ± 0.07 


0.38 ± 0.07 


0.5 ± 0.1 


GM Aur 


3.17 ± 0.20 


(20) 


3.29 ± 0.20 


1.29 ± 0.20 


1.6 ± 0.2 



nent of the surface density, and the average optical depth 
of the disk, respectively. 

For the surface-density profile we adopted a value of 
p = 1.5 throughout. The temperature profile index q is 
uniquely determined by the slope of the spectral energy 
distribution in the optically thick regime, given by the re- 
lation q = 2/(3 — «ir). We measured the spectral slope at 
far- infrared wavelengths (see Table [SJ, and derived tem- 
perature indices in the range of q « 0.5 — 0.7. Since A de- 
pends only weakly (logarithmically) on the average optical 
depth f , an order-of-magnitude estimate for this quantity 
will suffice for our purposes, f is a function of the opacity, 
mass, outer radius, a nd orientation of the disk (Eq. 16 in 
iBeckwith et al.llT990h . We estimate its approximate value 
by taking characte ristic numbers fo r these parameters, i. e. 
disk masses from iBeckwith et al. | dl990l) and disk radii 
from Table 2, for a circumstellar disk seen at random 
orientation. We found that f is of the order of ~10~ 2 . 
After inserting all parameters in the correction formula 
we derive A w 0.2. The optical-depth correction slightly 
increases the opacity indices. The final values are given 
in Table H 

4. Discussion 

One can use the millimetre spectral slope of circumstellar 
disks around pre-main-sequence stars to gain information 
on the characteristic dust grain size, provided the observed 
emission is largely optically thin. The magnitude of the 
dust opacity index indicates whether dust aggregates are 
small or large compared to the observing wavelength. The 
values are a robust measure of the characteristic grain 
size in the emitting region. At 7 mm, we largely observe 
the thermal emission of cold dust in the outer parts of the 
disk midplane. 

For dust grains with sizes of the order of a few 
tenths of a micron (<C A b s /27r), as present in the in- 



terstellar medium and in protostellar cores, the opac- 
ity index has been found to be ~ 2 llHildebrandl H983T: 
iDraine k Ledll984l:IOss"enkopf fc Hennindll994|) . For very 
large bodies (3> A b s /27r) that block radiation by virtue of 
their geometrical cross-section, the opacity is frequency- 
independent (grey opacity), = 0. 

Particles of about the same size as the observing wave- 
length of A = 7 mm, i. e. pebble-sized p articles, are ex- 
pecte d to be in an intermediate regime ijBeckwith et alJ 
l200(t) . Strictly speaking, one has to include the refractive 
index m = n + \k of the absorbing material. At millime- 
tre wavelengths, the optical consta nts for silicate grains 
are n ~ 3 and k ~ 0, thus l w I ~ 3 llLaor fc prainelll993t 
iHenning fc Mutschkd Il997t iMutschke et alJll998|) . Thus 
opacity indices <; 1 suggest the presence of dust parti- 
cles with sizes ;> A/2-7r|m|, i. e. with millimetre dimensions 
(and possibly even larger). 

After correcting for free-free radiation and optically 
thick emission, wc determined dust opacity indices in 
the range 0.5-1.6 for 10 sources where the circumstellar 
disk could be spatially resolved (Figure 0J. Six objects 
have value < 1, a robust indication of agglomerated dust 
particles in the millimetre size regime. Four other sources 
have values between 1.3 and 1.6, lower than the dust 
opacity index ~ 2 of submicron-sized particles as found 
in the interstellar medium. values smaller than 2 may 
also be explained by grain properties like shape, compo- 
sition, co nductivity, porosity, a nd crystallinity of the dust 
particles l|Henning et alJll995|) . Size, however, is the most 
important parameter influencing dust opacity indices. We 
deem unlikely the possibility that extreme dust composi- 
tions or particle structures cause th e observed low values 
llKoerner et al]ll 99,4 lDramfll2005h . 

We underline the importance of proper free- free correc- 
tion of the millimetre slopes, without which the opacity 
indices would be systematically smaller. The detection of 
emission at centrimetric wavelengths remains the only safe 
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method to estimate the contamination of free-free radia- 
tion to t he millimetre emission observed l)Testi et al .120011 
iNatta et al.ll2004 - 

5. Conclusion 

We presented continuum observations at 7 mm of a large 
sample of low-mass pre-main-sequence stars located in the 
Taurus- Auriga star-forming region with ages between ap- 
proximately 100 kyr and 3 Myr. Circumstellar material 
could be spatially resolved for 10 sources, allowing us to 
derive the size of their circumstellar disk. 

The deconvolved disk sizes (100-200 AU) effectively 
rule out the possibility that the 7 mm emission originates 
from small, optically thick disks containing submicron- 
sized dust particles. Knowing that the millimetre emis- 
sion is mostly optically thin, we could directly translate 
the spectral indices a mm measured in the millimetre part 
of the spectral energy distribution into the frequency de- 
pendence of the dust opacity. 

The power-law index (3 is indicative for the character- 
istic size of the radiating dust particles. From our 7 mm 
observation and other millimetre fluxes compiled from the 
literature, we found (3 values covering the range 0.5-1.6. 
More than half of the objects where the disk surrounding 
the star could be spatially resolved show clear evidence 
for the presence of millimetre-sized dust particles, demon- 
strating the operation of dust grain growth processes in 
disks around T Tauri stars. The remaining sources may 
not require large dust aggregates and are consistent with 
submicron-sized dust as found in the interstellar medium. 

Together with similar results for TW Hya, CQ Tau, 
and two Herbig Ae stars (HP 34282, HP 163296 ) 
l|Calvet et al.ll2002t iTesti et al.ll2003t INatta et al J 12004 - 
a picture of grain growth to pebble-sized particles in 
disks around low- to intermediate-mass pre-main-sequence 
stars is emerging. Further investigations of grain growth 
with existing and upcoming millimetre interferometers are 
needed to study where and how fast the building blocks 
of planetesimals are formed. 

Millimetre observations trace the cold outer disk mid- 
plane where most of the disk mass is locked up. The grain- 
size indicator (3 therefore only conveys information on that 
disk region. The surface layer of the inner disk is acces- 
sible through mid-infrared spectroscopy. The comparison 
of dust-size indicators from millimetre intcrfcromctry and 
mid-infrared spectroscopy, sampling different regions of 
circumstellar disks around T Tauri stars, may yield valu- 
able insights into the physical processes that are thought 
to control the formation of planetesimals and the evolu- 
tion of the circumstellar disk. Establishing the character- 
istic dust grain sizes at different locations in circumstellar 
disks holds the key to a better understanding of planet 
formation and disk evolution. 
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Fig. 4. Distribution of dust opacity indices. ISM-sized dust grains have (3 ~ 2, while particles much larger than the 
observing wavelength A=7 mm have a frequency-independent opacity law ((3 = 0). The populated intermediate regime 
is indicative for the presence of millimetre-sized dust aggregates. 



Table 5. Millimetre measurements and IRAS photometry compiled from the literature. The columns 9-q list the HPBW; note that the IRAS PSF sizes are 
given in arcmin. 



A 


Flux 


Error 


8b 


Ref. 


Flux 


Error 


9a 


Ref. 


Flux 


Error 


6b 


Ref. 


Flux 


Error 


6b 


Ref. 


Flux 
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Ref. 


(mm) 


(Jy) 


(Jy) 


(") 




(Jy) 


(Jy) 


(") 




(Jy) 


(Jy) 
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(Jy) 


(Jy) 


(") 




(Jy) 


(Jy) 
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RY Tau 






FT Tau 






DG Tau B 






DG Tau 






HL Tau 




0.012 


17.74 


0.027 


1' 


(a) 


0.36 


0.045 


1' 


(a) 


















9.74 


0.031 


1' 


(a) 


0.025 


26.48 


0.054 


1' 


(a) 


0.57 


0.037 


1' 


(a) 










30.39 


0.065 


1' 


(a) 


31.18 


0.073 


1' 


(a) 


0.060 


18.91 


0.066 


2' 


(a) 


0.82 


0.041 


2' 


(a) 










31.28 


0.039 


2' 


(a) 


76.26 


0.810 


2' 


(a) 


0.100 


13.50 


2.501 


4' 


(a) 


















39.86 
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4' 
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77.95 
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4' 
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0.036 


22 


(b) 
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22 
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22 
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0.07 
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18 
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12 


(g) 
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0.042 


12 


(g) 


0.961 


0.096 


12 
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0.229 


0.017 


11 


(c) 


0.130 


0.014 


11 
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11 
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0.879 


0.019 


11 
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2.0 
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5 


(d) 


















0.14 


0.03 


27.5 


(f) 


0.161 


0.017 


5 
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2.6 
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0.0017 


4 


(h) 


















2.7 










0.025 


0.0022 


3 


(e) 


0.0494 


r\ Ann O 

0.0088 


0.7 
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0.0020 


3 


(e) 
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0.6 


(i) 


3.0 
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0.0040 


6 
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0.0040 


3 


(j) 
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0.031 


1' 
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0.97 


0.034 


1' 


(a) 


0.78 
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1' 


(a) 










0.25 


0.031 


1' 


(a) 


0.025 


1.76 


0.043 


1' 


(a) 


1.32 


0.045 


1' 


(a) 


1.30 


0.046 


1' 


(a) 


4.07 


0.041 


1' 


(a) 


1.07 


0.042 


1' 


(a) 


0.060 


2.37 


0.069 


2' 


(a) 


1.39 


0.084 


2' 


(a) 


2.15 
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2' 


(a) 


6.33 
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2' 


(a) 


3.08 
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4' 
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0.039 


22 
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22 
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0.050 


22 
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0.032 


22 


(b) 
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0.03 


18.5 


(f) 
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0.02 


18.5 


(f) 


0.38 
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12 
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12 
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0.028 
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